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Gold-capped Janus particles immersed in a near-critical binary mixture can be propelled using il¬ 
lumination. We employ a non-isothermal diffuse interface approach to investigate the self-propulsion 
mechanism of a single colloid. We attribute the motion to body forces at the edges of a micronsized 
droplet that nucleates around the particle. Thus, the often-used concept of a surface velocity cannot 
account for the self-propulsion. The particle’s swimming velocity is related to the droplet shape and 
size, which is determined by a so-called critical isotherm. Two distinct swimming regimes exist, 
depending on whether the droplet partially or completely covers the particle. Interestingly, the 
dependence of the swimming velocity on temperature is non-monotonic in both regimes. 


The study of self-propelling synthetic colloids is an area 
of intense active research [1, 2]. The out-of-equilibrium 
directed motion of these colloidal microswimmers is 
maintained by a constant energy input which originates 
from their own activity. The directed swimming, coupled 
to the particle’s rotational diffusion, leads to a significant 
increase in the effective diffusion coefficient [3-5] and to 
complex collective behavior, such as dynamical phase- 
separation [6-8] and clustering [9-11]. Optimization of 
the microswimmers design is essential for realizing appli¬ 
cations such as targeted cargo and drug delivery, parallel 
assembly and scavenging of contaminants [1, 12, 13]. 

The design of synthetic swimmers requires an un¬ 
derstanding of the underlying mechanisms for the 
self-propulsion, e.g. self-diffusiophoresis [3, 14-18], 
self-induced electrophoretic flow, [19, 20] and self¬ 
thermophoresis [4, 21]. In many realizations, the particle 
motion is attributed to a microscopically thin boundary 
layer adjacent to the solid-fluid interface, which interacts 
with a self-generated held, such as electrical potential, so¬ 
lute concentration and temperature. Body forces within 
this layer give rise to an apparent slip velocity at the 
surface [22] while the huid outside the interfacial layer is 
considered force-free. Thus, the particle motion is com¬ 
pletely determined by the slip velocity distribution on 
its surface [17, 23]. However, this simple picture breaks 
down when the self-generated held extends to a region 
with a size similar to that of the particle. In this Letter 
we explore such a scenario of self-dihusiophoresis due to 
a local solvent demixing, leading to a complex swimming 
behavior arising from the coupling of the self-generated 
chemical potential gradients and the huid motion. 

We focus on a recently realized new class of swimmer 
consisting of Janus colloids immersed in a near-critical bi¬ 
nary mixture. Local heating of the colloid surface and the 
ensuing solvent demixing propels these particles, which 
exhibit fascinating individual and collective behaviour 
[5, 6, 24-26]. A similar system was studied by Araki 
and Fukai [27] but in their simulations heating is peri¬ 
odically applied to the whole mixture. In this work we 
study the self-propulsion mechanism of a locally heated 


Janus swimmer (illustrated in Fig. 1) and show that it 
is fundamentally different from motion driven by interfa¬ 
cial velocities. Rather, we hnd that the particle motion 
is linked to the how at the edges of a stationary demixed 
droplet, also far from the particle surface, and that it 
depends strongly on the solvent and particle properties. 
These advective effects are ignored in Wiirger’s recent 
study of a similar system [28]. 

We consider a micronsized spherical colloid with radius 
R immersed in an unbounded homogeneous near-critical 
binary mixture. Half of the particle is gold-capped and 
continuously heated by irradiation. The gold layer thick¬ 
ness is typically of the order of lOnm, much smaller than 
the colloid radius, and its thermal conductivity is much 
larger than that of the liquid and colloid. Thus, the cap 
forms an isotherm [21] with a temperature Tq > T^o, 
where T^o is the fluid’s ambient temperature far from 
the colloid. Local demixing of the fluid occurs if the tem¬ 
perature around the colloid increases above the (lower) 
critical temperature Tc of the mixture, of which both 
components are assumed to have a molecular volume 
with a ~ 3A. Consequently, temperature gradients in the 
demixed phase lead to gradients of the mixture chemical 
potential which give rise to a local body-force. From az¬ 
imuthal symmetry, no net torque will act on the colloid 
but since the temperature field is not spherically sym¬ 
metric, the demixed droplet exerts a net force on the 
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FIG. 1. Schematic illustration of a Janus particle immersed 
in a near-critical binary mixture at a temperature Too < 27 . 
Illumination results in To > Tc at the gold-capped hemisphere 
(red line). At steady state, a demixed droplet is nucleated 
around the translating colloid, with a fluid velocity Uoo at 
infinity in the frame-of-reference co-moving with the particle. 
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colloid in the axial direction. Hence, we set a cylindri¬ 
cal coordinate system in a frame-of-reference co-moving 
with a colloid placed at its origin. Our aim is to find 
the axial velocity of the fluid far away from the colloid: 
Uoo = —UooZ such that the colloid is force-free. The 
main result of this paper is that the self-propulsion is a 
result of forces perpendicular to the colloid surface, and 
cannot be attributed to an effective surface velocity par¬ 
allel to it, as is common in many scenarios [14, 17, 21, 22]. 

We investigate the steady-state during the ballistic mo¬ 
tion of the colloid at a time scale much shorter than its ro¬ 
tational diffusion time [24]. To this end we employ a well- 
established diffuse interface based approach [29], that 
couples the mixture order parameter ip G [—1/2,1/2], the 
mixture chemical potential ksTp, the fluid velocity UooV, 
and the scaled temperature 0 = {T — Too)/ {Tq — Too) by 


the dimensionless equations: 

dip/dt =—V ■ {Teipipv — Vp) , (1) 

M = -eVV + /'M, (2) 

V • V = 0 , (3) 

V-r = Vp-f Ch”^Ca"VVM , (4) 

her^dQ/dt = -V • (PeT0t’ - V0) . (5) 


Here, all lengths are scaled by R and time is scaled by 
R'^/D, where D is the mixture inter-diffusion constant. 

Eq. (1) is the convective Cahn-Hilliard equation gov¬ 
erning the composition dynamics, where the relative 
magnitude of the composition advective flux, ipv, and 
diffusive flux, DVp, is measured by the the composition 
Peclet number, Pe,^ = UooR/D. The chemical poten¬ 
tial Eq. (2) is obtained via the bulk free energy a^f = 
(x — 2)ip^ -I- 4ip^/3, where x ~ ^/T is the Flory inter¬ 
action parameter, supplemented by the square-gradient 
term which accounts for interfacial tension [30]. The 
interfacial width is characterized by e = x^h^) where 
Ch = a/R is the Cahn number. Eqs. (3) and (4) are the 
Stokes equations governing the fluid flow. Here, p and 
T = Wv + Wv'^ are the dimensionless fluid pressure and 
viscous stress tensor, respectively, scaled by rjfUoo/R, 
where p/ is the fluid viscosity. The last term in Eq. (4) 
is the body force due to gradients in the chemical poten¬ 
tial, it is proportional to the inverse capillary number, 
Ca = o^rjfUoo/kBT, which measures the relative mag¬ 
nitude of viscous and surface tension forces. Eq. (5) is 
the energy equation for the fluid. Here, Pe^ = UooR/ce 
is the thermal Peclet number, where a/ = kf/{pfCf) is 
the fluid thermal diffusivity, and pf, Cf and kf are the 
fluid density, heat capacity and thermal conductivity, re¬ 
spectively. The Lewis number, Le = af/D, is the ratio 
of thermal to mass diffusivity. 

The force, F, exerted on the colloid by the fluid is 
estimated by applying the divergence theorem to Eq. (4): 

F = 2tt J dc c [pi -I- H — t] • n , (6) 


where c = cos 'd, i9 is the polar angle and n is the outward 
unit vector normal to the surface. By recasting ipVp = 
V • H one obtains [29] the Korteweg stress tensor, H = 
[Tf'iT) - f - - epV^p] 1 -b eVipVip. 

The physical properties are taken to best mimic the 
experimental setup in Ref. [24]. Therefore, we use prop¬ 
erties of silica glass for the colloid and those of a crit¬ 
ical water-2,6-lutidine for the fluid. Since Tq — Too is 
only of the order of IK we can assume that all material 
properties are constant [30], see the caption of Fig. 2 
for the values used in the calculation. The assumption 
that D is constant is a somewhat crude approximation, 
since it vanishes as a power law close to Tc [31]. In the 
temperature window we examine D is of the order of 
10“^^ — 10“^^m^/s [32]. We therefore use an “effective” 
value of H as a parameter and examine its influence on 
the swimming. In addition, the exact value of the con¬ 
tact angle of each hemisphere in water-2,6-lutidine is un¬ 
known. We therefore use an indicative contact angle of 
9i = tt/A for the hydrophilic gold cap throughout this 
work such that p > 0 corresponds to a water-rich phase. 
We will explore the influence of 02 at the other hemi¬ 
sphere on the swimming. 

Taking R = 0.5pm, H = 4 x 10“^^ m^/s and a large 
Uoo = 20pm/s leads to Pe,^ « 0.2, Per « 10“^, Le « 10^, 
Ch ~ 10“^ and Ca « 10“^. Even for this extremal Uoo, 
PeT ^ 1 and we can safely neglect the heat advection in 
Eq. (5). In addition, because heat diffuses much faster 
than mass, Le 3> 1, the temperature adjusts almost in¬ 
stantly to a composition perturbation and we may ne¬ 
glect also the time dependence in Eq. (5), which leads 
to the heat equation V^0 = 0. The same argument 
holds for the solid where we also solve the Laplace equa¬ 
tion. Hence, the temperature distribution only depends 
on the thermal conductivity contrast of the solid and 
fluid, ks/kf, where ks is the colloid thermal conductiv¬ 
ity. Therefore, we solve the heat equation once, and use 
the resulting temperature distribution as an input for the 
relaxation of the velocity and composition. This relax¬ 
ation is dominated by the body force ipVp in Eq. (4) 
since Ca ^ 1. 

The colloid is placed at the origin of a large cylindri¬ 
cal domain of height 2 G [—L, L] and radius r G [0, L] 
with L = 500 [30]. We use symmetry boundary con¬ 
ditions (BCs) at r = 0. At the other edges of the do¬ 
main we impose, using the appropriate BCs, a mixture 
of critical composition at a temperature Too (0 = 0) 
with a velocity Uoo [30]. On both colloid hemispheres 
we impose a no-slip BC for the fluid: Vi = 0, where 
i = 1,2 denotes the capped and uncapped hemisphere, 
respectively. At the capped hemisphere we set 0i = 1 
while for the uncapped hemisphere we a have continuity 
of the heat flux, n ■ kfVQout = n ■ kg'^Qin- The first 
BC for the composition at the colloid surface is no flux: 
n-Vp = 0. The colloid has two chemically distinct solid- 
liquid interfaces, for which we assume an excess surface 
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free energy F^, of the form a^F^jj/ksT = f jiipdAi, 
where ji measures the difference between the microscopic 
short-range interaction of the two solvent components 
and the solid. The wetting angles are then imposed 
using: n ■ = — tan {tt/ 2 — Oi) |(1 — nn) ■ where 

in this so-called geometric formulation of the wetting BC 
7 i = cos9i/y/2 [33]. This BC has proved useful in simula¬ 
tions of moving contact lines where it is known to result 
in an effective slip through the diffusive fluxes between 
the phases [33-35]. Thus, fluid motion due to the inter¬ 
action with the surface is actually resolved even though 
the no-slip BC is imposed. 

Modern day computational resources coupled to spe¬ 
cialized mesh generation allowed us to resolve the fluid 
spatial distribution over 6 orders of magnitude, from the 
microscopic scale set by the interfacial width, 0(a), up 
to the mesoscopic scale set by the domain size L [30]. 
Steady State. The resulting composition of the mixture 
around a force-free swimmer (F = 0) for several Too and 
fixed To = Tc-I- 0.5K is shown in Fig. 2 (a)-(d). The solid 
red line in each panel is the contour of the reduced critical 
temperature, 0c = (Tc — Too)/(To — Too), which can ac¬ 
count for many features of the swimming. Demixing only 
occurs inside the region bounded by the 0c isotherm, 
where T > Tc. In Fig. 2 (a) we also show several other 
contours of 0 < 0c (dashed lines). Demixing within 
these contours will occur for fixed Tq > Tc and increasing 
Too [see Fig. 2 (b)-(d)] or for fixed Too and increasing Tq. 
In both cases 0c decreases and thus the droplet grows. 

Fig. 2 also reveals that a single water-rich droplet is 
nucleated at the particle surface. Within the droplet the 
composition is inhomogeneous; ip is maximal at the sur¬ 
face and decays rather smoothly to the bulk value ip = 0 
because of the temperature gradients and the proxim¬ 
ity to Tc. The demixed region is clearly distinct from 
the bulk phase, as can be seen from the velocity vec¬ 
tors. Inside the demixed region the velocity is very small 
implying the droplet effectively moves together with the 
particle. Strikingly, we find that no significant slip oc¬ 
curs at the surface, in contrast to the results in Ref. [28]. 
The fluid weakly circulates inside the droplet [30] and the 
overall flow pattern is similar to that of Stokes flow past 
a viscous droplet [36]. 

We distinguish between two droplet shapes: (i) when 
0c < 0COV the droplet partially covers the particle as in 
Fig. 2 (a)-(c) and (ii) for 0c > 0cov complete covering 
occurs as in Fig. 2 (d). The covering temperature 0cov is 
closely related to the uncapped pole isotherm 0p = 0c, 
for which the demixed region should first encompass the 
particle. For the solid-fluid heat conductivity contrast 
ks/kf « 3.5 that we use, 0p « 0.7. However, a thin 
demixed region at the pole is energetically costly, and in 
fact our numerical solution gives 0cov ~ 0.66 somewhat 
smaller than Op and a discontinuity in the pole composi¬ 
tion at 0COV In comparison, for ks/kf = 1 the lower heat 
diffusivity in the solid leads to Op = 0.5 [21]. Hence, the 
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FIG. 2. Steady state composition ip (right) and contours of 
the scaled pressure p/200 (left) around a heated Janus parti¬ 
cle in the xOz plane. The temperature of the heated cap is 
To = Tc + 0.5K. The bulk mixture temperature approaches 
Tc from (a) to (d) with Tc — Too equal to (a) 2.75K, (b) 1.5K, 
(c) I.IK and (d) 0.5K. The red curve is the critical contour 
0c. Arrows are vectors of the velocity (right) and body force 
[Eq. (4)] (left). In (a) the dashed lines are the contours 
0 = 0.75,0.69,0.6,0.5. Here, D = 4 x 10“^^ m^s, Si = 7r/4 
and 62 = '7r/2. For the physical properties of the 2,6-lutidine- 
water mixture we take a = 3.4A, Pf = 987 kg/m^ pf = 2 
mPa s, Cf = 4.3 kj/(kg K) and kf = 0.39 W/(m K) [37-39]. 
For the silica particle we used ks = 1.38 W/(m K). 


droplet shape and therefore the swimming behavior are 
both quite sensitive to the conductivity contrast, which 
we thus identify as an interesting engineering parameter. 

To maintain the steady-state shape of the droplet as 
dictated by the contour 0c, the composition diffusive flux 
(X Vp balances the convective flux according to Eq. (1). 
Therefore, the body force oc tpVp, is primarily large at the 
droplet diffuse boundary where advection becomes signif¬ 
icant, see the vector maps in Fig. 2 (a)-(c). The body 
force may be large within the droplet when Uao is small 
and the droplet internal dynamics also becomes signifi¬ 
cant, see Fig. 2 (d). Notice that the pressure contours in 
Fig. 2 are approximately perpendicular to pS/p because 
within the droplet V • r is small, and thus Vp oc pV/r to 
first-order at steady-state. Moreover, at the front of the 
droplet (with respect to the particle direction of motion) 
the diffusive flux must balance an advective flux toward 
the droplet whereas at the rear of the droplet the ad¬ 
vective flux carries the mixture away from the droplet. 
Thus, also the body force at the droplet edge acts in op¬ 
posite directions relative to the fluid flow, resulting in 
the pressure distributions of Fig. 2 (a)-(c) exhibiting two 
regions: (i) a p > 0 region at the droplet rear and (ii) a 
region of p < 0 near the three phase contact line. 
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FIG. 3. Swimming velocity vs. reduced critical temperature 
curves for three quenching temperatures Tq — Tc- The droplet 
completely covers the colloid for 0^ < ©cov ~ 0.66 whereby 
the different curves collapse onto one another. For 0c > 0cov 
the curves are distinct with the swimming velocity increasing 
with quenching. The velocity of the swimmers in Fig. 2 in 
indicated with arrows and the inset shows the corresponding 
force profiles along the colloid contour. 

The anisotropy of the droplet shape produces an 
anisotropy also of the pressure and body force within 
the droplet, which is responsible for the particle motion. 
Thus, the self-propulsion is a result of forces perpendicu¬ 
lar to the surface, which are affected by the flow far from 
the surface, and cannot be mapped to an effective surface 
velocity. In Fig. 2 (a)-(c), the resulting force-free motion 
is in the direction of the uncapped hemisphere whereas 
the completely covered colloid in Fig. 2 (d) moves with 
the capped hemisphere on the front. 

Swimming Velocity. Fig. 3 shows the swimming velocity 
Uao as a function of 0c for three values of the quench¬ 
ing To ~ Tc, which increases with the laser intensity in 
experiments. Uao strongly depends on 0c and is of the 
order of 0.1 — 10/rm/s, in agreement with experiments. 
Two swimming regimes exist, for complete coverage at 
0c < 0COV the swimming is independent of the quench¬ 
ing as 0c completely determines the demixed state. In 
contrast, for 0c > 0cov the curves are distinct with Uao 
increasing with quenching. This is because here 0c meets 
the particle surface at an angle different from the contact 
angle, which leads to a competition between the demix¬ 
ing and surface energies close to the three phase contact 
line, with the balance shifted in favor of the demixing 
as quenching increases. Experiments are performed at a 
single constant Too, showing that the swimming velocity 
increases with the laser intensity [24], but our calcula¬ 
tions, which explore a large range of Too, indicate this is 
not always the case. 

Fig. 3 also reveals that Uao jumps at 0cov and is non¬ 
monotonic with 0c in both regimes. To understand this 
behaviour we plot in the inset of Fig. 3 the z compo¬ 
nent of the force [cf. Eq. (6)] along the particle contour. 


z ■ c{pl -I- n — r) • n, as function of c, where c < 0 (c > 0) 
corresponds the (un)capped hemisphere. The labels in 
the inset and the arrow labels in Eig. 3 correspond to 
the swimmers in the panels of Fig. 2. When 0c and 
Uao are large, as in Fig. 2 (a), the force profile has a 
positive maximum in each hemisphere. This is because 
the droplet is small and hence a large body force exists 
close to the surface at both hemispheres. When 0c de¬ 
creases and the droplet grows, as in Fig. 2 (b)-(c), the 
positive maximum at the capped hemisphere becomes a 
negative minimum while the maximum at the uncapped 
hemisphere grows in magnitude and moves to a larger c. 
This is a result of the anisotropic shape of the droplet. 
While the 0c contour is distanced from the capped hemi¬ 
sphere it remains close to the uncapped hemisphere and 
also covers a larger portion of it. Thus, the body force lo¬ 
calized at the droplet edge becomes more significant near 
the three phase contact line rather than near the particle 
rear, thereby accounting for the minimum in Uao- For 
0c < 0COV, Uao first jumps to a larger positive value but 
decreases with 0c and eventually becomes negative, since 
now the droplet edge at the front is also distanced from 
the particle. Finally, Uao vanishes when 0c is small and 
the large droplet becomes quasi-spherical (Too —> Tc). 

The nucleation of a second, water-poor droplet, at the 
uncapped hemisphere is unfavored unless the hemisphere 
is very hydrophobic. The small volume of the demixed 
region around the uncapped hemisphere compared to the 
capped counterpart entails a relatively larger surface en¬ 
ergy penalty. Nonetheless, we find that for 62 ^ O.Stt a 
water-poor droplet does nucleate at the uncapped hemi¬ 
sphere. This is accompanied by a reversal of the swim¬ 
ming direction [30]. Assuming that the uncapped silica is 
hydrophilic [24, 40], this result is consistent with the ex¬ 
perimental observation that particles with a hydrophilic 
gold cap swim with their cap at the rear whereas for a 
hydrophobic cap it is at the front [24]. 

We also examined the dependence of Uao on the diffu¬ 
sion constant. For a swimmer with the same parameters 
as in Fig. 2 (b) but with D 10-fold larger or smaller, the 
resulting velocity is 111/rm/s and 0.36/rm/s, respectively. 
As expected, Uao increases with D because the diffusive 
current is able to compensate for larger advection while 
maintaining the droplet shape. Further work is required 
to understand the more realistic scenario where D may 
vary in space due to thermal gradients near Tc. 

In conclusion, we have shown that a locally heated 
Janus particle in a near-critical binary mixture is pro¬ 
pelled by the chemical potential gradients at the diffuse 
interface of a nucleated droplet, arising from the bal¬ 
ance of diffusive and advective fluxes. Therefore, the 
self-propulsion cannot be described by an effective sur¬ 
face velocity. We hope that our results will stimulate 
further experiments to uncover the details of the swim¬ 
ming mechanisms of these intriguing particles and possi¬ 
bly explore other microswimmers propelled by non-local 
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